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information processor 



(N 

o 

-t— > 
O 

O 



Guanru Peng 1 ' 2,3 , Guilu Long 1 ^ 2 and Raymond Laflamme 3 ' 4 
1 Department of Physics, Tsinghua University, Beijing 100084, China 
2 State Key Laboratory of Low- dimensional Quantum Physics, Tsinghua University, Beijing 100084, China 
3 Lnstitute for Quantum Computing and Department of Physics, 
University of Waterloo, Waterloo, Ontario, N2L 3G1, Canada 
4 Perimeter Institute for Theoretical Physics, Waterloo, Ontario, N2J 2W9, Canada 

Anyons have exotic statistical properties, fractional statistics, differing from Bosons and Fermions. 
They can be created as excitations of some Hamiltonian models. Here we present an experimental 
demonstration of anyonic fractional statistics by simulating a version of the Kitaev spin lattice model 
proposed by Han et al[f|] using an NMR quantum information processor. We use a 7-qubit system to 
prepare a 6-qubit pseudopure state to implement the ground state preparation and realize anyonic 
manipulations, including creation, braiding and anyon fusion. A ^ x 2 phase difference between the 
states with and without anyon braiding, which is equivalent to two successive particle exchanges, 
is observed. This is different from the -n x 2 and 2ir x 2 phases for Fermions and Bosons after two 
successive particle exchanges, and is consistent with the fractional statistics of anyons. 

PACS numbers: 03.67.Lx,05.30.Pr,03.67.Pp 
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I. INTRODUCTION 

In 3-dimensional space, indistinguishable particles 
obey Fermi-Dirac statistics (Fermions), or Bose-Einstein 
statistics (Bosons). For both Fermions and Bosons, upon 
the exchange of two indistinguishable particles, the sys- 
tem wave function gains a n or 2n phase change. How- 
ever, when restricted to 2-dimensional space, particles 
appear to obey fractional statistics @- This means that 
when two indistinguishable particles in 2-dimensional 
space are exchanged, the system wave function gains 
a statistical phase change, ranging continuously from 
to 27r. Those quasiparticles are defined as anyons. 
Anyons can be grouped in abelian and nonabelian 
anyons. Abelian anyons are particles that realize 1- 
dimensional representations of braid groups. In nature, 
abelian anyons are believed to exist and be responsi- 
ble for the fractional quantum hall effect (FQH)0-Q. 
Nonabelian anyons are particles that behave as multi- 
dimensional representations of braid groups. They are 
critical in topological quantum computing, for example, 
in the Kitaev fault-tolerant quantum computation mod- 
els @, Recently, the interest in anyons is enhanced 
by the developing field of quantum computing because of 
their potential ability to impliment fault-tolerant quan- 
tum computing architecture. 

Several theoretical schemes have been proposed to di- 
rectly observe the fractional statistics associated with the 
anyon braiding motion 0, Gil- These schemes are 

mainly grouped into two approaches: the first is pro- 
posed to be realized in FQH systems, and the second 
makes use of the Kitaev models. In FQH systems, it 
is difficult to directly observe anyonic fractional statis- 
tics, and to introduce or resolve individual anyons [l6j 
when compared with the schemes using the Kitaev spin 
lattice models. Experimental demonstrations in photon 
systems using the Kitaev spin lattice models have been 



realized [171 Ho]. However the anyons are not protected 
from local noise and there is no explicit particle inter- 
pretation of the excitations [l5j because the background 
Hamiltonian vanishes in such photon systems. In con- 
trast, the background Hamiltonian can be simulated in 
the nuclear magnetic resonance (NMR) systems. 

In our work an NMR quantum information processor is 
used to demonstrate the anyon braiding scheme proposed 
by Han et al.[l[ in the smallest Kitaev system utilizing 
6 qubits. The six-body ground state preparation, anyon 
excitations, and anyonic braiding operations are realized 
using a 7-qubit molecle in liquid state NMR. By compar- 
ing the two final states, of which one is obtained after the 
anyon creation, braiding, and fusion processes while the 
other doesn't undergo such processes, the phase differ- 
ence, which is mapped into a frequency change of NMR 
spectrum peaks in our experiment, can be observed. 



II. THE KITAEV k x k SQUARE LATTICE 
MODEL 

The first Kitaev spin lattice model@ is a k x k square 
lattice on the torus, containing qubits on each of the 
bonds (Figure Q} (here we define the bonds as the mini- 
mal lines forming the lattice) . The total number of qubits 
is 2fc 2 . The spin lattice contains vertices and faces. A ver- 
tex v is the intersection of four bonds. A face / means 
the surface with boundary defined by four bonds. We 
can then define an Hamiltonian as 



where 
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A v (Bf) represents the 4-body interactions belonging to 
the qubits which live on the vertex v (or the face /). For 
all the vertices and faces, the ground state \^ ground) for 
the Hamiltonian Hk satisfies 

A v \^ ground) — \ Aground) , Bf | ^ ground) — \ ^ ground) • (3) 

The ground states are 4-fold degenerate. They form a 
protected subspace G 

G = {|0 e N,A v \0 = |0,% 10 = 10, for all v and /}. 

(4) 

N is the Hilbert space of the 2k 2 qubits. This is the 
definition of the toric code, which is a special kind of 
stablizer code@. A v and Bf are its stabilizer operators. 
Because of the periodic boundary conditions, for each 
qubit j, (Tj and crj appear twice in different A v 's and 
B/ s. It can be easily obtained that 

n v A v = l,U f B f = 1. (5) 

If a state does not satisfy several (for example, n) of the 
A v \^)} = j*) and Bf\<&) = |>3>) constraints, it's an excited 
state with n elementary excitations (or quasiparticles) . 
The relationships in Eq. [5] implies that the quasiparticles 
should appear in pairs. 







— •— 






4 


► < 


► 4 
• 


> \ 














—•— 










• 








• 




• 





FIG. 1: Illustration of the first Kitaev model. The first Kitaev 
spin lattice model is a k x k square lattice on the torus. There 
is a qubit on each bond. The operaters A v and Bf act on the 
qubits of the vertex v and face /, respectively. 



qubit ee 




(a) (b) 



FIG. 2: Illustration of anyon creation and braiding opera- 
tions, (a) The creation of two m (e) anyons by the operation 
Xmm (Zee)- Xmm means ax operation to qubit mm, and Zee 
means a z operation to qubit ee. The two m anyons are local- 
ized on face fl and face /2. The two e anyons are localized 
on vertex vl and vertex v2. (b) The braiding motion of an 
m anyon around an e anyon by operations X11X12X13X14. 
Qubits 11, 12, 13, and Z4 are the qubits along the braiding 
path. 



The states with quasiparticles (excitations) are excited 
states. (Figure [Ha)) 

Since two m (e) particles at the same site annihilate, 
the m (e) particle can be moved by applying a x (a z ) 
operations along the path (Figure [2](b)) . A braiding op- 
eration is to move an m (e) particle around an e (m) 
particle along a closed circle path, which is equivalent to 
two successive particle exchanges [19]. For Fermions and 
Bosons, states do not change after two successive particle 
exchanges. For the m and e particles, it has been shown 
that after a braiding operation, the global state gains a 
I x 2 phase change 16], which is different from Fermions 
and Bosons. Therefore the m and e particles are anyons 
that obey fractional statistics. 



If a o x operation is applied to a qubit, for example, 
qubit mm, in a ground state, the state wave function 
is I*) mm = X mm \^ ground) {X m m means the a x opera- 
tion to qubit mm). It satisfies Bf\\ i i>) mm = — \^) mm 
and Bf 2 \f) mm = -|*)mm- /I and /2 are the two 
faces next to qubit mm. That means two quasiparti- 
cles have been created at those particular locations. The 
two quasiparticles can be considered as two "defects" lo- 
calized on faces fl and /2. They are called m parti- 
cles. Instead, if a cr z operation is applied to a qubit, 
for example qubit ee, in a ground state, the state wave 
function is |*) ee = Z ee \§> grouml ) (Z ee means the a z op- 
eration to qubit ee). It satisfies A v \\§>) ee = — |\E') ee and 
A V 2\^)ee = — \^)ee- This also means two quasiparticles 
occur. Here, vl and v2 are the two neighboring vertices 
which are connected by the bond qubit ee lives on. The 
two quasiparticles can also be considered as "defects" lo- 
calized on vertices vl and v2. They are called e particles. 



III. THE SIX-QUBIT KITAEV SPIN LATTICE 
MODEL AND THE EXPERIMENTAL SCHEME 

The minimum amount of qubits needed to implement 
the smallest version of the periodic Kitaev model for 
anyon braiding operations is eight. However by aban- 
doning periodic condition, the spin lattice model can be 
extended from a square lattice to any planar graph [20jj . 
and an anyonic model can be found with six qubits where 
we can demonstrate braiding statistics [lj. The graphic 
structure of the six-qubit model is shown in Figure [3](a). 
The Hamiltonian of the system is 



H 6 = -A x -A 2 -B 1 -B 2 -B 3 ~ B 4 , (6) 
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The ground state of the six-qubit Kitaev spin lattice is 

Aground) = |(|000000> + 1111000) + 1110111) + 1001111)). 

(7) 

Because the boundary conditions have changed, the 
ground state is not degenerate anymore. The ground 
state can be created from a six-qubit graph state shown 
in Figure [3[b) [l[ . A graph state is a type of multi-qubit 
state represented by a graph with the vertex set V and 
the edge set E, and is defined as 

\G) = n {i ^ E u id \+)® v , (8) 

where the operater U^j is the controlled-cr z operation be- 
tween the qubits i and j, and |+) = ^(|0) + |1)). The 
six-qubit graph state corresponding to the graph in Fig- 
ure [3Jb) is 

\G e ) = i:,J ] :;i:; A J:J-,<. + + + + + +). (9) 
The ground state of the six-qubit Kitaev spin lattice is 

^ground) = 0\G 6 ), (10) 

where O = IHHHHI, I is the identity operater, and 
H is the Hadamard operater. The ^ground) can be pre- 
pared by first preparing a six-qubit graph state \Ge) , then 
implementing an O operation. This gives a method to 
prepare the ground state. 

For the ground state of the six-qubit system in Figure 
l^a), if a a x operation is applied to qubit 4, a pair of m 
particles are created on its two neighboring faces, and if 
a a z is applied to qubit 3, a pair of e particles are created 
on its two neighboring vertices. By applying successive 
a x operations on qubits 6, 5, 3, and 4, one m particle 
moves in a loop around one e particle. After such a 
braiding operation, the global wave function will obtain 
a phase factor of -1. By observing this phase change, one 
can verify the fractional statistics of anyons. 

Han et al.[l| give the basic circuit for ground state 
preparation, anyon creation, anyon braiding, and anyon 
fusion. It should be noted that if one does the braid- 
ing operation to a state with a pair of e particles and a 
pair of m particles (Figure Hfb)), the phase change af- 
ter braiding is a global phase, which cannot be observed 
in experiments directly. In the scheme proposed by Han 
et al.Q, the anyon creation step is realize by a x and 
•J<7z = e l ~e~ % ~° z instead of a x and<7 2 . The a x operation 




(a) (b) 

FIG. 3: (a) The six-qubit Kitaev model and its braiding loop. 
A pair of m (e) anyons are created by the operation X4 (Z3) 
and the braiding operation is realized by X6X5X3X4. (b) The 
graph state that is equivalent under local unitary operations 
to the ground state of the Hamiltonian in (a). The corre- 
sponding qubits in the NMR spin system are also labeled at 
each vertex. 

creates a pair of m particles. The ^fal operation creates 
a superposition between the states with and without a 
pair of e particles. Therefore after the anyon creation 
step, the state of the system is 

l*) = ^(lV'i) + W), (11) 

where is a state with a pair of m particles only, 
\ip2) is a state with a pair of m particles together with a 
pair of e particles. will not change after a braiding 
procedure, because no e particles exist. \ip2) will obtain 
a phase factor of -1 after a braiding procedure. Therefore 
the total wave fuction becomes 

|*') = i(|^)-|V2». (12) 

In this way, the phase change caused by braiding opera- 
tion becomes a local phase factor in front of \tji2) an d is 
observable in experiments. 

The fusion operation is realized by applying ^JaZ and 
a x . In such case, if indeed there is a statistical phase 
change | x 2 acquired after braiding, the state after fusion 
is \^ ground)- This means the e particle pair and the m 
particle pair are both fused. 

In our experiment scheme, yfolT 1 and a x operations 
are performed as the fusion step (Figure [SJ . With the 
statistical phase change | x 2 introduced by a braid- 
ing operation (two successive exchanges between anyons) , 
the state after the fusion step should be \^ excited) — 
o ground) ■ Without the statistical phase change § x 2, 
the state after the fusion step should be | $ gr0 und) ■ There- 
fore, by observing the difference between the state after 
the fusion step and the ground state, we can demonstrate 
the fractional statistics of anyons. 

IV. EXPERIMENTAL IMPLEMENTATION 

In the experiment, 13 C-labeled trans-crotonic acid dis- 
solved in d6 acetone was used. The system contains 7 



4 



qubits. The molecule and its parameters are described 
in Figure |U 

To implement the anyon braiding, we first prepared 
the molecule in the labeled pseudopure state with a devi- 
ation matrix of the form pi — Zh\^c\^m^C2^c^H2^cz 
where = |0)(0| and Z is the Pauli matrix a z . Wc 
chose HI as the label qubit, and CI, M, C2, C4, H2, C3 
as qubits 1, 2, 3, 4, 5 and 6, to match to Figure (3Jb), 
using the neighboring couplings for shortening the gate 
operations in implementation. 
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FIG. 4: Characteristics of the molecule of trans-crotonic 
acid[2lj, 22]. The chemical shifts (diagonal elements) and J- 
coupling constants (off-diagonal elements) are given in Hz. 
The spin-lattice and spin-spin relaxation times Tl and T2 
are listed at the bottom. The chemical shifts are given 
with respect to reference frequencies of 700.13 MHz (hydro- 
gens) and 176.05 MHz (carbons). The three hydrogen nu- 
clei in the methyl form a spin-3/2 group M. After using a 
gradient-based subspace selection, this group acts in its spin- 
1/2 subspace 21]. Therefore, M can be used as a single qubit. 
We thus have 7 qubits, including the two protons and four 
carbons. 



FIG. 6: The quantum network for the experiment without 
anyonic manipulation. It serves as a comparison with the 
experiment in Figure [S] 



from which one can observe the phase change from the 
spectrum of the label qubit. These two experiments il- 
lustrated in Figure [5] and Figure [5] were carried out for 
comparison. 

The wave fuctions of the labeled states (a), (b), (c), 
(d), (e) in Figure[Sl and (f), (g) in Figure[B]are as follows: 

l*a) = ^(|0ci0 M 0c 2 0c40//20c3) + \IciImIc2®C$H2®Cz) 

+ |1ci1m0c21c41h21(73) + |0ci0Mlc2lc4.1ff2lc3/) 

— ground) 



\^b) = -(|0ci0M0c2lc40ff20c3) + « I lei 1 M lc2 ^CA^H2^Ci) 



+ |1ci1m0c20c41// 2 1 

1 
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Xc?4(-y=(e^ Aground) 
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:<*))) 
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FIG. 5: The quantum network for the experiment with any- 
onic manipulation. H represents the Hadamard operation. 
Z = a z , X = a x . F90 is the read pulse, F90 = e~ i2 i" y . 



— 2 (|0ci0Mlc20c4lff 2 lc3/ + «|lcilM0c20c4lff 2 lc3) 
+ |1ci1m1c 2 1-C40h20c3) + «|0ci0A/0c 2 lc40ff20c3/) 



I* 



excited 



))) 



X ci {-j={e 1 ^ I Aground) 



2 (|0ci0A/0c20c40_ff 2 0c3) — |lcilAfl-C20c40/f20c3 

|1ci1m0c21c41j?21c3) — |0ci0m1-C21c41.(Y21-C3/) 
iZc2\^ ground) = * | ^excited) 



Figure [S] and Figure [5] show the circuits for the ex- 
periments with and without anyonic manipulation, re- 
spectively. There is a part indicated as "Measurement" 



I*, 



iV2 



(|0ci0A/lc20c40ff 2 0c3) + |1cj,1m1c21c41h21c3)) 
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I*/) = 2(I°ci0m0c20c40 H2 0c3) + |lcilMlc20c40 ff2 0c3; 

+ \^C1^mOc2^C4^H2^C3) + \^Cl^M^-C2^-CA^H2^-Cz)) 
— I ^ ground) 

-^-(|0ci0jl/0c20c40if20c3) + I Icil M^C2^Ci^ H2^C3j) 

(13) 

In the experiment with anyonic manipulation, \^ a ) is 
obtained after the ground state preparation. It is the 
ground state of the six-qubit Kitaev model. After anyon 
creation, the wave function of the state is \^b)- After 
anyon braiding, is transformed to \^ c ). \^b) and 
l^c) are the superpositions of \4>i) and \if)2)- is a 

state with a pair of m particles. 1^2) is a state with a 
pair of e particles and a pair of m particles. Due to the 
e anyons present in \ip2), it has a phase change acquired 
when the m particle braiding around the e particle. This 
causes the difference between \^fb) and \^f c )- After anyon 
annihilation, the state goes to the excited state \^fd)- In 
order to observe the phase change in relatively simple 
spectra directly, the state \^d) is transformed to |\& e )- 
For the state l^e), by observing the label qubit, a 2-peak 
spectrum can be obtained. The left peak in the spectrum 
corresponds to 1 111111) while the right peak corresponds 
to 1001000). 

In the experiment without anyonic manipulation, the 
ground state preparation and the measurement proce- 
dures are the same as the experiment with anyonic ma- 
nipulation. The state |^/) equals l^a)- The state \^ g ) 
is almost the same as l^e), except that the qubit C2 is 
|0) in \ty g ) (the global phase is ignored). For the state 
l^g), by observing the label qubit, a 2-peak spectrum 



can be obtained. The left peak in the spectrum corre- 
sponds to 1 110111) while the right peak corresponds to 
1 000000). It is straightforward to conclude that the dif- 
ference is caused by anyon braiding, illustrated in Figure 
[5] Because m and e anyons do not obey integral statis- 
tics, after two successive exchangings they obtain a phase 
factor, which is mapped into a frequency change of peaks 
in our experiment. 

At the end of the measurement part, the states of 
qubits HI and C2 are exchanged via a swap gate. There- 
fore, C2 becomes the label qubit. The final results are 
obtained by implementing a ir/2 read pulse to C2. It 
should be mentioned that the J-coupling constant be- 
tween C2 and H2 \Jc2,H2\ — 0.66Hz is the smallest of 
the couplings between C2 and the other nuclei. | Jc2,H2\ 
can be resolved in our experiments, which means all 
the 64 peaks of C2 spectrum are resolved (Figure [7]). 
\Jc2Mi\ — 155A2Hz is the largest J-coupling constant 
of C2. The braiding induced changes in frequencies of 
peaks should equal | Jc2,m|, which can be easily and ex- 
plicitly observed in the experiments. 

The gates used in the ground state preparation were 
realized by combining single-qubit rotations and evo- 
lutions of the J-coupling constants between the neigh- 
boring qubits, while all the anyonic manipulations were 
realized by single qubit rotations [U, [25|. The single- 
qubit rotation pulses were generated using the GRAPE 
algorithm[26j for HI and H2, and were standard Isech- 
shaped r.f. pulses for M and C1-C4. The J-coupling evo- 
lutions were realized by implementing refocusing pulses. 
We combined all the pulses using a custom-built software 
compiler, which numerically optimizes refocusing pulses 
and minimizes the errors due to undesired J-coupling 
evolutions [13, The duration of the pulse sequences 
shown in Figure[B]and Figure[S]was 195.1ms and 250.2ms, 
respectively. 



V. EXPERIMENTAL RESULTS 

The final results for the experiments are shown in Fig- 
ure [71 Figure 0(a) shows the superposed spectra between 
the simulated thermal state spectrum of C2 and the ex- 
perimental pseudopure state spectrum. It shows the ex- 
perimentally realized 1 000000) peak (peak o). It should 
be mentioned that there is an antiphase peak (peak w) in 
Figure[7]Ja). This antiphase peak was caused by the label 
qubit HI, which was in the Z state. Figure EKb) displays 

the spectrum of the state 1 ^(|0ci0M0m0c40,f/ 2 0c3) + 
|lcilM0mlc4l#2lc3)), and Figure Efc) displays the 
spectrum of the state 

1 2^(|0ci0a/1hi0c40//20c3) + |lcilj\/lmlc4lff2lc3))> 
both generated by observing C2, whose deviation density 
matrix was X = <j x . Peaks i (-18.3i/z), j (-137.2iJz), s 



(137. 1Hz) and t (18.2Hz) correspond to states 1 110111), 
1000000), 1111111) and |001000), respectively. The sum 
intensity of the two dominant peaks, both in Figure [TJ^b) 
and Figure Etc), is about 0.7, normalized using the in- 
tensity of peak o. 

There is a 155 AH z — \Jc2.hi\ shift between peaks i 
and s, j and t (Figure [7]). This frequency difference be- 
tween the peaks in the two spectra was caused by the 
process of anyonic manipulation, demonstrating that af- 
ter the braiding operation, the state with e and m anyons 
acquired a phase change S = (§ + rj) * 2. Here 2rj is the 
deviation of the experimental phase change from | x 2. 

We suppose the experimentally realized labeled state 
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FIG. 7: (a) The superposed spectra of the theoretical C2 thermal state spectrum (red line) and the experimental pseudopure 
state spectrum (black line). There are 64 peaks in the thermal state spectrum, each corresponding to a computational basis 
state, (b) The experimental spectrum corresponding to the experiment in Figure HJ] It has two dominant peaks, i and j. (c) 
The experimental spectrum corresponding to the experiment in Figure [5] It has two dominant peaks, s and t. The amplitude 
of peak o in the experimental pseudopure state spectrum in (a) is taken as reference to normalize the experimental signals 
shown in (b) and (c). On the right are the zoomed-in spectra for the peaks o, i, j, t, and s. The states which the experimental 
peaks correspond to are labeled on top of each zoomed-in spectrum. There is a Jc2,h\ distance between peaks i and s, j and 
t. p, q, u, and v are the small peaks that have the same frequencies as those of peaks s, t, i, and j, respectively. 



(a) in Figure [5] and labeled state (f) in Figure [6] were 



I*/') = |V) = a\^ ground ) + ^exczted) error)- 

(14) 

The above expression implies the ground state prepa- 
ration was not perfect. (3\^ excited) was responsible for 
peaks p and q. 'y\^f e rror) was responsible for the peaks 
other than peaks i, j, p and q. (Figure [7(b)) The experi- 
mentally realized labeled state (g) in Figure [6] was 

IV) = ^-[«(|0ci0m0c 2 0c40h20 C 3) 

+ \^C1^M^C2^C4Xh2^-G3)) 

+/3(|0ci0 m 1c20c40h 2 0c3) 

error / ■ 

(15) 

| — | can be determined from the peak intensities (denoted 
as r) of the experimental spectrum shown in Figure[jjb). 



,0. 



T p + T « =0 .i8 ±0,09. 

r 4 + r, 



(16) 



The experimentally realized labeled states (d) and (e) 
in Figure [5] were 

|*/) = ie"' [(-a sin 77 - /3 cos ri)\^ ground ) 

+ (a cos X] - /3 sin 77) | ^ excited) } + 7 1 <A error / 

= ie %V [a' \9 ground) + 13' \^ excited)} + l^error) ■ 

(17) 



IV) = 



y/2ie ir > 



;a'(|0ci0 M 0c20c40ff 2 



C3 



+ |lcilM0c 2 lc4lff2l-C3)) 
+ /3'(|0ci0A f lc20c40tf 2 0c3) 
+ |lcilAflc2lc4lff2lc3))] " 



(18) 



Here 77 - s 



a 



-a sin 77 — (j cos 77, ft 1 = a cos 77 



/3 sin 77. a' \ty g roun d) was responsible for peaks u and v. 
7 1 terror) was transformed from error) via anyonic ma- 
nipulation. j\4>error) was responsible for the peaks other 
than s, t, u and v. (Figure [7(c)) 

I ^7 1 can be determined from the experimental spec- 
trum shown in Figure [7(c) . 



a', /r u + r„ 



= 0.24 ±0.06. 



(19) 



1 /?" v r s + r ( 

Combining Equation [16] and Equation I19[ we can obtain 



I — I - 1-1 

tan n = 2_ = 0.O6 ± 0.03. (20) 

l + l|l*lv#l 



77 = 0.06 ± 0.03 = (0.02 ± 0.01)tt and the phase change 
S = (f + 77) * 2 = (0.52 ± 0.01)tt x 2. This agrees with 
the prediction of the fractional statistics. 
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VI. DISCUSSION 

The signal loss mainly came from the spin-spin relax- 
ation and pulse imperfection. Through comparing the 
signal intensities of the simulation with and without the 
T2 effects, we estimate that T2 effects contributed to 
about 15% of the loss of signal. Comparing the experi- 
mental results (peak intensity ~0.7) with the results of 
the simulation with T2 effects (peak intensity ~1), we 
estimate that imperfections of the implementation of rf 
pulses caused an additional approximate 15% signal loss. 

The values of r\ in the simulation with and without T2 
effects are 0.07 and 0.06, respectively, which match well 
with the experimental data (77 = 0.06 ±0.03). The differ- 
ence between the two r\ values in the simulation is small, 
which means dccohercncc did not contribute much to the 
deviation of 5 from | x 2. The deviation was mainly 
caused by imperfections of the refocusing protocols and 
the implementation of rf pulses. 



VII. CONCLUSION 

In summary, we demonstrated the anyonic fractional 
statistics, using a 7-qubit NMR system. This is the 



first demonstration on topological quantum computing 
using nuclear spins. One advantage of our experimental 
scheme is that we can use the same technique to simulate 
a 9-qubit Kitaev spin model mentioned by Han et al.[l[ 
using an NMR system with more spins, so that we can 
do a demonstration that anyonic operation is robust to 
different braiding paths and thus taking an important 
step towards showing the fault tolerance properties of 
the Kitaev model. 
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